Determining neutron-capture cross sections via the surrogate reaction technique 
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Indirect methods play an important role in the determination of nuclear reaction cross sections 
that are hard to measure directly. In this paper we investigate the feasibility of using the so-called 
surrogate method to extract neutron-capture cross sections for low energy compound- nuclear 
reactions in spherical and near-spherical nuclei. We present the surrogate method and develop a 
statistical nuclear-reaction simulation to explore different approaches to utilize surrogate reaction 
data. We assess the success of each approach by comparing the extracted cross sections with a 
predetermined benchmark. In particular, we employ regional systematics of nuclear properties in 
the 34 < Z < 46 region to calculate (n, 7) cross sections for a series of Zr isotopes, and to simulate 
a surrogate experiment and the extraction of the desired cross section. We identify one particular 
approach that may provide very useful estimates of the cross section, and we discuss some of the 
limitations of the method. General recommendations for future (surrogate) experiments are also 
given. 



PACS numbers: 24.10.-i, 24.60.Dr, 25.40.Lw, 98.80.Ft 



I. INTRODUCTION 

Nuclear reaction cross sections are often difficult to 
measure directly. This is particularly true for reactions 
relevant to applications in nuclear astrophysics, since ra- 
dioactive nuclei play an influential role in many cosmic 
phenomena but cannot be easily studied in the labora- 
tory. Information on these nuclei, and on their relevant 
cross sections, is needed to improve our understanding of 
the processes that shape our universe. While indirect 
methods for determining direct-reaction cross sections 
have become very popular in recent years [l[ , compound- 
nuclear reaction cross sections are typically determined 
purely theoretically. 

For nuclei with mass A > 30, statistical-reaction model 
(Hauser-Feshbach) calculations are widely used to esti- 
mate cross sections that have not been measured. These 
calculations require input data such as masses, onc- 
and two-particle separation energies, properties of res- 
onances, level densities, optical potentials for particle 
transmission coefficients, and gamma strength functions. 
This data should be constrained by measurements wher- 
ever possible, but for the thousands of nuclei involved 
in astrophysical environments one has to rely on global 
phenomenology or, alternatively, on microscopic nuclear 
theories. 

In this paper, we explore the possibility of using an 
indirect method, the surrogate nuclear reaction method, 
for obtaining (n, 7) compound-nuclear reaction cross sec- 
tions. A simplified version of the method, which com- 
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bines experiment with reaction theory to obtain cross 
sections for reactions that proceed through a compound 
nucleus, was first used in the 1970s [2, 9 to extract (n, /) 
cross sections for various actinides from transfer reactions 
with t and 3 Hc projectiles on neighboring nuclei, followed 
by fission. A modern version of the approach was used 
by Petit et al. to study the 233 Pa(n, /) reaction cross sec- 
tion using a (^He,p) transfer reaction [J. More recently, 
Burke et al. Q and Plettner et al. [|| constructed ratios of 
decay data from surrogate experiments on two different 
uranium isotopes, and used that information to extract 
the 237 U(n, /) cross section. The surrogate ratio method, 
as applied to actinide nuclei, was examined in much de- 
tail by Escher and Dietrich [7| . These efforts have shown 
that the surrogate method is a very useful tool for pre- 
dicting various (n, /) cross sections in the actinide region, 
and there is no a priori reason why the method should be 
limited to studies of fission. In principle, the surrogate 
method can be applied to any reaction that proceeds via a 
well-defined, equilibrated compound state; but its great- 
est potential value lies in applications that involve unsta- 
ble isotopes. Among the unanswered mysteries about the 
nature and evolution of our universe is the origin of the 
heavy elements. Much effort is currently being devoted 
to exploring possible paths for the nucleosynthesis (such 
as the s and r processes) and astrophysical environments 
that can produce the elements between iron and uranium. 
Of particular interest in the context of the s process are 
(n, 7) reactions on branch point nuclei, unstable isotopes 
with a life time long enough to allow the reaction path 
to proceed by either neutron capture or (3 decay. In prin- 
ciple, these isotopes are ideally suited for investigations 
using the surrogate method since they are located next 
to stable elements that can be used as targets in the 
surrogate experiment. Surrogate approaches for other 
neutron-induced reactions have been considered as well. 
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Early experiments [8[ were carried out to assess the feasi- 
bility of using the surrogate technique to determine cross 
sections for (n, a) and (n,p) reactions on nuclei in the 
mass-90 region. These experiments highlighted several 
issues that needed to be addressed in order to extract 
reliable cross sections from surrogate measurements. 

The experiments mentioned above were analyzed un- 
der the simplifying assumption that the decay probabil- 
ities are independent of the particular spins and pari- 
ties of the compound-nuclear states that are occupied in 
the neutron-induced as well as in the surrogate reaction. 
This assumption, which is known as the Wcisskopf-Ewing 
limit 0, [13] , is not always valid and its application needs 
further exploration. We will investigate the validity of 
the Weisskopf-Ewing approximation for (n, 7) reactions 
involving spherical, or near-spherical, targets in the Zr 
region. The fact that the spin-parity ( JII) distributions 
in the decaying compound nucleus can be very different 
in the n-induced and the surrogate reaction is referred to 
as the " JIT population mismatch" . For s-process branch 
points, e.g., low-energy neutrons bring in very little angu- 
lar momentum, while direct reactions leading to the same 
compound nucleus can produce very different angular- 
momentum distributions. This leads to challenges for 
extracting information from a surrogate experiment, as 
was already recognized from the early experiments in the 
mass-90 region [8(. 

Recently, Younes and Britt [111, [HJ demonstrated that 
taking into account the JII population mismatch can 
have a significant effect on the extracted results. They 
revisited the data from the original surrogate transfer- 
reaction induced fission measurements from the 1970s, 
and employed a simple direct-reaction model to account 
for the angular-momentum population difference between 
the neutron-induced and direct reactions. They were 
able to improve on earlier results for (n, /) cross sections 
for various Th, U, Np, Pu, and Am isotopes. For 235 U, 
Younes and Britt reproduced the known fission cross sec- 
tion for the J = 7/2 + ground state and were able to pre- 
dict the fission cross section for the isomeric l/2 + state 
at 77 eV, which to date has not been measured directly. 

In this paper we investigate the feasibility of applying 
the surrogate method to extract low-energy (n, 7) cross 
sections on mass ~ 90 — 100 nuclei. There are essentially 
two different sources of uncertainty inherent in the sur- 
rogate method: (i) Insufficient knowledge of the decay 
pattern for the relevant compound nucleus, which must 
be supplemented by reaction modeling; (ii) Insufficient 
knowledge of the spin-parity distribution of the decaying 
compound nucleus. We will present calculations that il- 
lustrate the effect of these two sources of uncertainty on 
cross sections extracted from surrogate experiments. We 
will consider and compare different strategies of utiliz- 
ing the surrogate data for obtaining unknown compound- 
nuclear reaction cross sections. 

We will present calculations for a range of zirco- 
nium isotopes and in particular wc perform surrogate 
experiment simulations to study the extraction of the 



91 Zr(n, 7) 92 Zr cross section. However, our discussions 
are more general and can in principle be applied to most 
spherical and near-spherical nuclei in the mass 90-100 
region. This mass region is not only of importance for 
understanding the s-process nucleosynthesis path, but 
also encompasses the majority of light fission fragment 
nuclei of the binary fission yield distribution. Conse- 
quently, there is great interest in obtaining cross sec- 
tions for neutron-induced reactions involving nuclei in 
this region. A particularly interesting application ap- 
pears in nuclear astrophysics: 93 Zr and 95 Zr are both 
s-process branch points and are believed to be produced 
in AGB stars. It has been proposed [l3[ that the rela- 
tive abundance of 94 Zr and 96 Zr, as measured in prcsolar 
grains [Til ], depends on the efficiency of mixing between 
the H- and He-burning shells of AGB stars. Thus, better 
knowledge of the neutron capture cross sections of 93 Zr 
and 95 Zr can lead to a diagnostic tool to probe stellar 
interior physics. 



II. METHOD OF THE STUDY 

In the present work we use a statistical-reaction model 
simulation to assess whether the surrogate method can 
be employed to extract low-energy (n, 7) cross sections. 
First, we demonstrate how the typical level of uncertainty 
in the Hauser-Feshbach parameters affects the cross sec- 
tions obtained in a purely theoretical approach. The 
cross section that best represents the available data will 
later on serve as a benchmark for the surrogate method. 
The combination of models and parameters used to pro- 
duce this benchmark cross section is referred to as the 
"Reference Decay Model" and will be considered as the 
most realistic description of the true decay of the com- 
pound nucleus. We simulate the impact of having insuffi- 
cient information about the "true" decay by considering 
variations of the decay model. We also use the same set of 
decay models when studying the decay of the compound 
nucleus populated through a surrogate reaction. This 
type of theoretical approach has a number of distinct 
benefits: (i) We are able to access physical quantities 
that are not directly observable in an experiment, such 
as spin-parity dependent branching ratios for individual 
exit channels, (ii) We can alter key quantities such as 
level densities and gamma strength functions and carry 
out sensitivity studies, (iii) Performing simulations will 
allow us to identify the main limitations of the method 
and to quantify the precision that can be expected. 

All calculations are carried out with a modified version 
of the Hauser-Feshbach code Stapre [ll| [l(| ■ In order 
to account for the fact that a direct reaction produces a 
jn distribution in the residual nucleus that is different 
from the one associated with the desired neutron-induced 
reaction, we have modified the code so that it is suitable 
for surrogate-reaction studies. In particular, wc included 
an option to allow the JII distribution of the first com- 
pound nucleus to be read in from a file rather than cal- 
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culated from entrance-channel transmission coefficients. 
It is therefore possible to specify an arbitrary JII distri- 
bution for a given compound nucleus and to predict the 
decay of the nucleus. In addition, we have implemented 
the specific models for level densities and photon strength 
functions that we are using for this work (see Sec. lIV Al) . 

The theoretical framework for the surrogate method is 
described in Sec. lIIII Statistical- reaction calculations are 
performed using recently developed regional systematics 
for nuclear properties in the 34 < Z < 46 region. These 
systematics, and the resulting estimated (n, 7) cross sec- 
tions for a range of zirconium isotopes, are presented in 
Sec. [IVJ In Sec. [V] we use the reference decay model 
obtained from the regional systematics to simulate a sur- 
rogate experiment. We discuss three different approaches 
for utilizing surrogate data to extract low-energy neutron 
capture cross sections. We perform sensitivity studies 
and discuss theoretical and experimental challenges for 
the surrogate method. Concluding remarks and recom- 
mendations are given in Sec. I VII 

III. STATISTICAL REACTION THEORY AND 
THE SURROGATE IDEA 

A. Hauser-Feshbach theory 

The formalism appropriate for describing compound- 
nucleus reactions is based on the statistical Hauser- 



wherc s x and l x are, respectively, the exit channel spin 
and the orbital angular momentum between the residual 
nucleus and the ejectile. Jq and lie denotes the spin 
and parity, respectively, of a state in the residual nucleus 
C. The excitation energy is given by Uq, and energy 
conservation gives Uc = U — S c — e x . The terms in the 
sums are furthermore restricted by parity conservation 
and spin-coupling conditions. The particle transmission 
coefficients, T x (e x , l x ), can in general depend on the com- 
plete set of quantum numbers that define a channel, but 
the dependence is here limited to energy and orbital an- 
gular momentum. For gamma-decay channels, the sum 
over l x , s x is replaced by a sum over electromagnetic mul- 
tipolcs, and the transmission coefficients accordingly re- 
fer to these multipoles. pc{Uc, Jc, lie) denotes the level 
density of the residual nucleus. The quantum numbers 
and quantities for all open channels, x', appearing in the 
denominator, are defined analogously. Finally, although 
the above formula does not indicate it in the interest of 



Fcshbach theory |f0|, |I7[ ■ In this section we summarize 
the Hauser-Feshbach theory in a form suitable for appli- 
cation to surrogate experiments that will be discussed in 
Sec. IIII Bl Let us consider the case of a reaction leading 
from an entrance channel a + A (denoted a) via an inter- 
mediate compound-nuclear state B* to an exit channel 
c+C (denoted x) ■ The energy-averaged cross section for 
this reaction, as a function of the center-of-mass energy 
e a in the incoming channel, is written as 

v« x (e a ) = a ^ J, n)G x (C7, J, IT), (1) 

where we have assumed that the Bohr hypothesis of inde- 
pendence between formation and decay holds separately 
for each value of total angular momentum and parity, 
JII, of the compound nucleus. The formation cross sec- 
tion, a a (e a , J, II), for particular values of JII in the com- 
pound nucleus, as well as the branching ratio of decay 
into channel x, G X (U, J, II), are both energy-averaged 
quantities. The excitation energy U of the compound 
nucleus is related to the separation energy S a of particle 
a by U = S a + e a . 

For radiative capture reactions at low energies one is 
usually interested in the cross section for producing a 
particular isotope. Therefore, we will define G X (U, J, II) 
as the branching ratio of decay into channel \ integrated 
over all bound states of the residual nucleus 



(2) 



simplicity, Stapre actually treats each step of the de- 
cay completely and correctly; in particular, for primary 
gammas populating levels that are still above the par- 
ticle thresholds, Stapre calculates the competition be- 
tween particle and gamma emission in dealing with fur- 
ther steps in the cascade. 

The assumption of full independence between forma- 
tion and decay of the compound nucleus can be re- 
laxed by the introduction of width fluctuation corrections 
w ax (e a , J, II) . These aim to correct for correlations be- 
tween the widths of the incoming and outgoing channels. 
In general, the width fluctuations will only be important 
for very small energies, where the effect will enhance the 
elastic channel and simultaneously deplete all other chan- 
nels due to flux conservation. We will demonstrate in 
Sec. HVI that. for the cases of interest here, the depiction 
of the (n, 7) cross section is < 10% below I McV, and 
negligible above that. 



GJU,J,IL) 



Sz„,s Y S/J c ,n c Jo 



u-s c 



dU c T x (e x ,l x , )pc{U c , Jc, n c ) 



E x < Ei x „ Sx , Ej c „n c „ Io~ Sc ' dUc>T x ,(e x ,,l x ,)p c >(Uc>, Jc^c) 

I 
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B. The surrogate method 

Here we outline how the cross section of Eq. ([1]) can be 
determined by a combination of theory and experiment 
in the surrogate method. In a surrogate experiment, the 
compound-nuclear state B* is produced via an alterna- 
tive ("surrogate") direct reaction d+D — > b+B*, and the 
decay of B* into the desired channel x is observed in coin- 
cidence with the outgoing particle b. In the following we 
assume that the ejectile, b, is observed at a particular an- 
gle flp and for simplicity we suppress the angular depen- 
dence of the observed quantities. The energy-differential 
cross section for the direct reaction, constituting the first 
step of the surrogate reaction sequence, is then 



das 
dep 



das p , 



(3) 



where ep is the center-of-mass energy in the outgo- 
ing channel b + B* (denoted /?), and JTI on the right- 
hand side denote different spin-parity combinations of the 
states of the compound nucleus populated via the incom- 
ing channel 5 of the surrogate reaction. This decompo- 
sition must be provided by a reaction-model calculation. 
Making the non-trivial assumption that B* damps into a 
fully equilibrated state, we can define the probability for 
producing the compound nucleus at excitation energy U 
with spin and parity JII 



^(£/3, J,n) 



(4) 



For a fixed ejectile angle the reaction kinematics will 
give a straightforward relation between the projectile and 
ejectile energies, and the excitation energy U of the com- 
pound nucleus. 

In a surrogate experiment the ejectile b is measured 
in coincidence with an appropriate observable that tags 
the channel \ of the desired reaction (e.g., a fission frag- 
ment for neutron-induced fission reactions, or an emitted 
gamma-ray for capture reactions). The experimental ob- 
servable of interest is therefore the probability of this 
coincidence 



p x (u) = E F w( u > J > n ) G x C 17 . J < n ) > 



(5) 



J. II 



where we have combined the formation probability of 
Eq. (J3J) with the branching ratio of Eq. @ . As mentioned 
above, the production of equilibrated compound-nuclear 
states following the direct reaction is a non-trivial issue 
that requires further attention. For the purpose of the 
present study we assume that such states have been pro- 
duced and that these states subsequently decay according 
to the branching ratio G (U, J, II). 

In principle, width fluctuation corrections should be 
incorporated in the above expression. However, if the 
formation of the compound nucleus, represented by 



F S a(U, J, II), is entirely independent of its decay, i.e. has 
none or negligible contribution to the total decay width, 
then the width fluctuation corrections can be applied to 
the branching ratios alone. Surrogate reactions such as 
inelastic scattering, are examples of this type of reac- 
tion. In this work, we will neglect width fluctuations in 
the modeling of the surrogate decay probabilities, but we 
will introduce them in the final step where wc extract the 
desired (n, 7) cross section, see Eq. (JH)). 



C. The angular momentum mismatch 

The relationship between Eqs. (TT]) and (O constitutes 
the Hauser-Feshbach formulation of the surrogate reac- 
tion method. While the standard Hauser-Feshbach for- 
mula expresses a cross section in terms of products of 
formation cross sections and decay branching ratios, the 
experimental observable in the surrogate approach is a 
coincidence decay probability which is expressed in terms 
of formation probabilities and the decay branching ra- 
tios. Although these two expressions contain the same 
branching ratios, G X (U, J,H), they are weighted differ- 
ently because the JTI distributions are different in the 
two reactions. In fact, this mismatch can be quite signif- 
icant. As we will show in Sec. El the use of the surrogate 
technique to extract low-energy (n, 7) cross sections is 
particularly challenging since the neutron will bring in 
very little angular momentum to the compound nucleus 
compared with that typically brought in by the D(d, b)B* 
direct reaction. A theoretical challenge in this case is 
therefore to determine the Jli distribution F S g(U, J, IT) 
so that the branching ratios G X (U, J, TV) can be extracted 
from the observed coincidence decay probability, P (U), 
sec Eq. ([5]). If this decomposition is well determined, the 
branching ratios can be inserted in Eq. (TTJ) together with 
a formation cross section a a (e a , J, n) calculated with an 
optical potential, to yield the desired a ax {e a ) cross sec- 
tion. At this final stage we can also introduce width 
fluctuation corrections, represented by the factor w ax in 
the following expression. This procedure for analyzing a 
surrogate experiment can be outlined in diagram form 



measured 



calculated 



extracted 



P X (U) =E Fsp(U,J,U) GJU,J,U) 



j. 11 



deduced 



E ^a,J,n) G x (c/ 5 j,n) u; QX (£ a ,j,n) 



calculated 



calculated 



(6) 

At this point we can remind ourselves that we have as- 
sumed that the ejectile, b, is observed at a particular 
angle Clp, or alternatively that the experimental data is 
integrated over a limited range of angles, and that we 
have neglected all angular dependence in the formulas 
above. However, the formation probabilities F sg (U, J, H) 
are most certainly angle-dependent, and it is therefore 
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desirable for surrogate experiments to be carried out over 
a wide range of ejectile angles in the D(d, b)B* reaction. 
Obtaining the same a ax for the various angles provides 
an important consistency check on the procedure. 

An important simplification occurs when the branch- 
ing ratios of Eq. do not depend on J and EL This re- 
moves most of the model dependencies from the surrogate 
analysis since the angular-momentum mismatch becomes 
irrelevant. This limit is known as the Weisskopf-Ewing 
limit Q . A number of conditions must be satisfied in or- 
der for this limit to be applicable 0, SEBl- For example, 
the energy of the compound nucleus must be sufficiently 
high so that all channels into which it can decay are dom- 
inated by integrals over level densities. This condition is 
not satisfied for low-energy neutron radiative capture and 
the dramatic simplification of the surrogate analysis that 
occurs in the Weisskopf-Ewing limit, see e.g. Ref. [H, 
cannot be utilized. The breakdown of the validity of the 
Weisskopf-Ewing approximation at low incoming neutron 
energies is illustrated in the next section, see in particu- 
lar Fig. |4j Furthermore, the spin-dependence of the level 
densities in the relevant channels has to be of the sim- 
ple form p(U, J) oc (2J + 1). This condition is satisfied 
if the spin J is smaller than the spin-cutoff parameter 
CT cu t in the relevant level-density formula. However, it 
is known that the Weisskopf-Ewing limit is still a useful 
approximation at somewhat higher spins. 



IV. NEUTRON CAPTURE ON ZR ISOTOPES 

In this section we present calculated (n, 7) cross sec- 
tions for a range of Zr isotopes, obtained by applying 
Hauser-Feshbach theory with decay models adjusted to 
regional systematics. The purpose is threefold: Firstly, 
we want to illustrate the quality of the developed sys- 
tematics and present the uncertainties connected to cer- 
tain input parameters. Secondly, these calculations pro- 
vide a prediction for the unmeasured 95 Zr(n, 7) 96 Zr cross 
section which is of importance for stellar modeling ef- 
forts, see e.g. Refs. [H, [H, EH- And finally, the cal- 
culation presented here for 91 Zr(n, 7) 92 Zr will serve as a 
reference when discussing the feasibility of using surro- 
gate experiments to deduce (n, 7) cross sections. We use 
our statistical-reaction decay model to generate gamma- 
decay probabilities and we are able to perform sensitivity 
tests by varying the most relevant input parameters. 



A. Decay model and regional fit 

Neutron capture cross sections modeled with the 
Hauser-Feshbach formula are most sensitive to the pho- 
ton transmission coefficients and nuclear level densities. 
For the modeling of capture reactions with stable tar- 
gets, these two inputs can usually be constrained by 
experimental data. In particular, the average total s- 
wave radiation width, (r 7 ) , may be used to determine 



the overall normalization factor associated with the El 
photon strength function while measured neutron s-wave 
resonance spacings, Dq, can be used to fix level-density 
parameters at the neutron separation energy. For cases 
where these quantities have not been measured one may 
attempt to describe them by systematics. 

A recent effort has been made to develop regional sys- 
tematics (as opposed to global prescriptions) for these 
two quantities for nuclei spanning the range 34 < Z < 46 
(selenium through palladium) [20| . The goal in develop- 
ing these systematics was to model several neutron- and 
charged-particle reactions on isotopes of yttrium, zirco- 
nium, niobium, and molybdenum. Here we make use of 
the findings of this work. Below we summarize the es- 
sential features of the model. More details can be found 
in Ref. [H. 

The level density is described by a backshifted Fermi 
gas combined with a constant temperature form at low 
excitation energies [2l], ■ It is assumed that the parity 
distribution of nuclear states is equal, i.e. for a given 
parity n, p(U, J, n) = p(U,J)/2. The spin-dependence 
of the level-density is described by 



f(U,J) 



2J + 1 
~2^~ 2 



■ exp 



-(J 



2a 2 



(7) 



(8) 



where the spin-cutoff parameter is given by 

a 2 = 0.01496AA 5 / 3 ^. 

In our analysis we fix A = I . The level-density parameter 
a, introduced above and also appearing in the Fermi gas 
formula, is given an energy dependence 



a{U, Z, N) = a(Z, N) 



1 + SW{Z,N) 



1 - exp [-7(17 - A)] 



U-A 



(9) 



where a(Z, N) is the asymptotic value for large U, 
5W(Z, N) is the shell correction as defined in [22| . 
(U — A) is the backshifted energy, and 7 is a constant 
factor. Values for the backshifts are taken as the aver- 
age difference in binding energies of neighboring nuclei, 
as described in (23[. Also taken from [23[ are the val- 
ues a(Z, N) = 0.1337(Z + N) - 0.06571(Z + N) 2 ' 3 and 
7 = 0.04884 MeV -1 . For nuclei with measured D n , a 
shell correction 6W(Z, N) is determined so that it will 
reproduce the measured level spacing. For other nuclei, 
the shell correction is based on the systematics shown 
in Fig. [TJ These systematics were determined by mak- 
ing two least-squares quadratic fits (one for either side of 
the closed N = 50 neutron shell) to the extracted shell 
corrections 



6W{Z, N) = cqN 2 + Cl N + c 2 , 



(10) 



where 



co = 
-0.114399 
-0.0421006 



ci = 
9.44901 
5.45665 



C2 = 

-188.821, 
-171.285, 



for N < 50 
for N > 50, 
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valid for 34 < Z < 46. The error bars on the extracted 
shell corrections reflect uncertainties in the measured D . 
The level-density parameters related to the constant tem- 
perature form used at low excitation energies are fixed 
by the known discrete spectrum and by requiring the 
two level density formulae to match tangentially at an 
energy U m . This matching energy may be adjusted to 
provide the best possible fit to the low-lying level struc- 
ture without affecting the Fermi gas portion of the level 
density. Such adjustments were made individually for 
each nucleus. 
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FIG. 1: (Color online) Systematics for the shell correction 
used in the Fermi gas description of the level density. See 
text for details. 

The systematics for the average total s-wave radiation 
width were determined by assuming a simple dependence 
on the mass and s-wave resonance spacing: 



<rx s = M- 



b] [c\og 10 (D ) + d] 



(ii) 



Making a least-squares fit to measured (r 7 ) (in meV) 
and Do (in keV) taken from the compilations of |24[ yields 
the coefficients a = -5.928, b = 343.8, c = -42.43, and 
d = 343.8. These systematics were then used to generate 
a list of suggested values of (r 7 ) for nuclei in the range 
34 < Z < 46. The list includes the available experimen- 
tal values. Systematic values were chosen based on mea- 
sured Dq where available, or on values of Dq calculated 
from the regional level density systematic shown in Fig.[TJ 
The line shape of the El photon strength function is de- 
scribed using a simplified version of the extended gener- 
alized Lorentzian (EGLO) model [2j|. The EGLO model 
corresponds to a Lorentzian with an energy-dependent 
width. We are using a simplified version in which the 
width is set to depend only on the nuclear temperature 
at the energy of the decaying state and not on the energy 
of the final state. The most important feature, however, 
is that the magnitude of the strength function is nor- 
malized to reproduce the average total s-wave radiation 
width. 

Since we are interested in modeling reactions for a 
relatively large region of isotopes we use the global 
nucleon-nucleus optical-model potential by Koning and 
Delaroche \2m for the calculations of the nucleoli trans- 



mission coefficients. It has been shown [26| that this pa- 
rameterization gives a very satisfactory fit to measured 
total cross section data and low-energy observables, such 
as the s-wave strength function, that are relevant for our 
application. 



B. Neutron capture cross sections for Zr isotopes 

Our calculated (n, 7) cross sections for the selected 
range of Zr isotopes 90-95 Zr are shown in Fig. [51 The 
experimental data for the quantities that constrained our 
decay model is summarized in Table [I] As noted earlier, 
the two most important ones are the average total s-wave 
radiation width, (r 7 ) , and the measured s-wave reso- 
nance spacing, D . For the two radioactive isotopes, 93 Zr 
and 95 Zr, for which no such experimental data exist, we 
have used values obtained from the regional systematics 
described in the previous section. The error bars given for 
these cases are estimates based on typical experimental 
uncertainties for other isotopes in this region. Error es- 
timates for our calculated (n, 7) cross sections are shown 
as dashed and dash-dotted lines in Fig. [5] They were ob- 
tained by repeating the calculations using the upper and 
lower limits of the s-wave radiation width and level spac- 
ing, respectively. Our calculated results are compared to 
available experimental data from the EXFOR/CSISRS 
database [23, [H, El, EL EL EL El, El, HI and we find a 
satisfactory agreement for the isotopes that are well stud- 
ied and are characterized by high level densities, which 
implies that the statistical-reaction treatment is well jus- 
tified. However, the resonance spacing in 91 Zr is rather 
large and signatures of individual peaks are observed in 
the 90 Zr(n, 7) 91 Zr experimental data which suggests that 
the statistical treatment might not be appropriate for this 
case. In TableUwe also present the Maxwellian-averaged 
(71,7) cross section (MACS) at k B T = 30 keV. Our re- 
sults are compared to values recommended by Z.Y. Bao et 
al. [Ill , which are based on an evaluation of available ex- 
perimental and theoretical results. 



V. THEORETICAL ANALYSIS OF DECAY 
DATA FROM A SURROGATE EXPERIMENT 



In this section we explore different approaches for uti- 
lizing decay probabilities measured in a surrogate exper- 
iment to extract low-energy neutron capture cross sec- 
tions. We perform our simulations for the 91 Zr(n, 7) 92 Zr 
reaction. Let us first take a closer look at some of the de- 
tails of the statistical-reaction calculation introduced in 
the previous section. The decay model for 92 Zr described 
in Sec. IIV Al will act as our reference model for these stud- 
ies, and the corresponding (n, 7) cross section serves as 
the reference (or "true") cross section that we seek to 
extract from our simulated surrogate experiments. 
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FIG. 2: (Color online) Calculated (n, 7) cross sections c omp ared to experimental data for six Zr isotopes: 90 Zr, 91 Zr, 92 Zr, 93 Zr, 
94 Zr, and 95 Zr. Experimental data are from Refs. H3, HI [H, HJ 00], H SJ H3, IH . Dashed and dash-dotted lines indicate 
error bars associated with uncertainties in the s-wave resonance spacing Do and the s-wave average radiative width (r 7 )o, 
respectively. The error bars for the 93 Zr and 95 Zr (n, 7) reactions are estimates based on typical experimental uncertainties for 
other isotopes in this region. 



TABLE I: Maxwellian-averaged (71,7) cross section at thermal energy ksT = 30 keV calculated for different Zr targets. The 
results from this work are compared to the recommended values by Z.Y. Bao et al. [3(|. Also shown in the table are: Q 
(Q-value for the reaction), No (number of known discrete states in the compound nucleus that form a complete spectrum), U m 
(matching energy of constant temperature and Fermi-gas level-density regions), Do (level spacing for s-wave resonances at the 
neutron separation energy), and (F 7 ) (s-wave average radiative width). Experimental data are from the RIPL-2 database [24| . 
whereas unmeasured quantities are obtained from our regional systematics. 





90 Zr(n, 7 ) 


91 Zr(n, 7 ) 


92 Zr(n, 7 ) 


93 Zr(n, 7 ) 


94 Zr(n, 7 ) 


95 Zr(n, 7 ) 


Q [MeV] 


7.194 


8.635 


6.734 


8.221 


6.462 


7.856 


N D 


41 


6 


3 


9 


1 


8 


U m [MeV] 


3.167 


4.658 


2.054 


3.058 


2.0 


2.695 


D [keV] 


6.0 ± 1.4 


0.55 ±0.10 


3.5 ±0.8 


0.160 ±0.015 


3.2 ±0.8 


0.26 ± 0.04° 


(r 7 ) [meV] 


130 ± 20 


140 ± 40 


135 ± 25 


164 ± 40° 


85 ±20 


144 ± 40" 




Maxwellian- 


averaged (n, 7) cross 


section (in n 


b) at thermal energy ksT 


= 30 keV 




This work " 


26 ±7 


62 ± 19 


41 ± 11 


198 ± 44 


27 ±9 


117 ±34 


Bao et al. [36] 


21 ± 2 


60 ±8 


33 ±4 


95 ± 10 


26 ± 1 


79 ± 12 c 



"There exists no direct experimental information for this observ- 
able. The given value is based on the systematics from our regional 
fit. 

''The error bar given here is \j Aj, + A 2 ,, where Ao r are the 

errors due to the uncertainties in Do and (r 7 )o, respectively. 
c This recommendation is based purely on theoretical estimates. 
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TABLE II: A summary of the decay models for 91 Zr(n, 7) 92 Zr 
that are used in this work to investigate the sensitivity to 
certain key quantities. 





<r 7 > 

[meV] 


D 

[keV] 


SW 
[MeV] 


MACS at 
30 keV [mb] 


Reference 


140 


0.55 


-1.743 


62 


Decay model 1 


280 


0.55 


-1.743 


115 


Decay model 2 


140 


0.45 


-1.237 


74 



A. The 91 Zr(n,7) 92 Zr reaction 

Neutron transmission coefficients calculated with the 
global optical-model potential by Koning and De- 
laroche [26[ are shown in Fig. E)Ja) for several partial 
waves on a 91 Zr target. For / > 0, these are the appropri- 
ately weighted averages of the coefficients with j = Z ± i . 
The rapid increase of the p-wave transmission coefficient 
with increasing neutron energy is characteristic of all iso- 
topes in this mass region, and corresponds to the well 
known 3p giant single-particle resonance. An interest- 
ing consequence of this feature is that low-energy neu- 
tron absorption (above approximately 10 keV) populates 
predominantly the JITstates of the compound nucleus 
that can be reached by either s-wavc or p-wave capture. 
The effect is clearly illustrated in Fig.[3Jb) where the JH 
population of 92 Zr, following the absorption of 30 keV 
neutrons on 91 Zr(5/2 + ), is shown to be dominated by 
1 _ ,2~,3 _ ,4~ and 2+,3+ states. 

It is important to realize that there is generally a very 
small number of open decay channels following the ab- 
sorption of low-energy neutrons on spherical, or near- 
spherical, targets. Other particle-decay thresholds usu- 
ally appear above the neutron separation energy, and in- 
elastic channels are closed when the first excited states of 
the target are above the incoming neutron energy. This 
leaves only two possible decay paths: compound-elastic 
scattering leading back to the target ground state, and 
gamma deexcitation of the compound nucleus. We now 
focus on the branching ratio defined in Eq. @. In this 
special case, the sum over open channels (%') appear- 
ing in the denominator contains only two terms. The 
important quantities arc therefore the neutron transmis- 
sion coefficients discussed above and the product of the 
gamma transmission coefficients with the level density 
of the compound nucleus. The latter quantity must be 
integrated over all states that can be reached by a pri- 
mary gamma transition. For low-energy neutrons, a sin- 
gle gamma deexcitation is usually enough to bring the 
compound nucleus below the neutron-decay threshold, 
which will inevitably lead to a continued gamma cascade 
down towards the ground state. As noted previously, 
the Stapre code correctly computes the gamma cascade 
even when neutron emission is possible within the cas- 
cade. 

Fig. 5] shows the gamma branching ratios for 92 Zr as 




FIG. 3: (Color online) (a) Energy- and ^-dependence of neu- 
tron transmission coefficients on a 91 Zr target, (b) JTl pop- 
ulation of the 92 Zr compound nucleus following neutron ab- 
sorption on 91 Zr(5/2+) at e„ = 30 keV. 



a function of JTI and neutron energy. For comparison 
we also plot the total gamma decay probability which 
we define as the 92 Zr production cross section divided 
by the neutron absorption cross section. These results 
are obtained with width-fluctuation corrections turned 
off. The results shown in this plot demonstrate a feature 
of nuclei that is particularly striking near closed shells, 
namely a dramatic dependence of the compound-nucleus 
gamma decay branching ratios on JIT at low energies 
above threshold. For some JIT values the probability 
for gamma decay is close to one, whereas for others it 
is on the order of 10 -3 . Clearly we are very far from 
the Weisskopf-Ewing limit and certain approximations 
that have previously been used in the analysis of surro- 
gate experiments cannot be used here. The reason for 
the dramatic dependence on JIT is the behavior of the 
neutron transmission coefficients, shown in Fig. [3Ja), in 
combination with the small number of discrete states at 
low excitation energies in the initial nucleus 91 Zr. Those 
compound- nuclear states in 92 Zr that can reach an en- 
ergetically allowed state in 91 Zr by the emission of an 
s- or p-wave neutron have a very small gamma-decay 
probability. For the other states, the branching ratios 
are determined by the competition between I > 2 neu- 
tron transmission coefficients and gamma transmission 
coefficients. It turns out that the gamma-decay chan- 
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nel usually dominates although the gamma transmission 
coefficients are still very small. It is also important 
to note the convergence of the curves at higher ener- 
gies E„ > 3 MeV. This demonstrates the onset of the 
Weisskopf-Ewing limit with the resulting approximate 
independence of JTI for the branching ratios. In the 
Sec. IV CI we return to this issue and discuss how we can 
utilize the observed behavior for our purposes. 

Before we discuss the sensitivity to properties of the 
nuclear decay model, we introduce the notation a^(e n ), 
a lbs( £ n), and G™ f {U, J,II). These quantities are the re- 
sults of our reference calculation and they denote, respec- 
tively, the radiative neutron-capture cross section, the 
neutron absorption cross section, and the gamma branch- 
ing ratios of Eq. fTJ). 

By introducing variations from the reference decay 
model we will be able to study the effects of our lack 
of knowledge of the true decay model, which is one of 
the two major uncertainties associated with the surro- 
gate method. As previously mentioned, the modeling of 
neutron capture reactions is most sensitive to the nu- 
clear level density at, and somewhat below, the neutron 
separation energy, as well as to the El photon strength 
function at small gamma energies (^2 — 3 MeV). In 
order to quantitatively investigate this sensitivity we in- 
troduce two decay models, "Decay model 1" and "Decay 
model 2", in addition to the "Reference Decay model" 
corresponding to the parameter set from the regional sys- 
tematics. The most important parameters for the three 
models are summarized in Table [TT1 together with the re- 
sulting MACS at 30 keV. The salient feature of Decay 
model 1 is a factor of two increase in the magnitude of 
the El photon strength function. Such uncertainty in the 
size of the gamma-transition strength is not unrealistic 
when moving to unstable isotopes for which no exper- 
imental data exists. The (n, 7) cross section is almost 
directly proportional to the magnitude of the El photon 
strength function, as is clearly seen in Fig.[5]where the ra- 
tio of the modeled capture cross section to the reference 
calculation is shown. The logarithmic energy axis was 
introduced in order to cover a large energy range while 
still focusing on the relevant small energies. For Decay 
model 2 we decreased the level spacing D n for s-wave 
resonances by modifying the shell correction SW in the 
formula for the level-density parameter a, Eq. We 
choose this particular way of modifying the level density 
since SW is the free parameter in our regional system- 
atic. The < 20% decrease of the level spacing at the 
separation energy implies a > 20% increase of the level 
density, and again we observe a proportional increase in 
the capture cross section. Finally, we present the results 
obtained when turning off the width-fluctuation correc- 
tions. As expected, turning off width-fluctuation correc- 
tions results in an increase in the capture cross section. 
From Fig. [S] we infer that it is approximately a 10% ef- 
fect at small energies and completely negligible at higher 
energies (> 2 MeV). 




I 1 I 1 I 1 I 1 I 1 I 

1 2 3 4 5 
e n [MeV] 

FIG. 4: (Color online) Gamma branching ratios as a function 
of spin, parity and excitation energy of the decaying state in 
92 Zr. The energy is given as the equivalent neutron energy 
e n = U — S n - The total gamma-decay probability following 
neutron absorption, as described in the text, is shown as a 
solid line in both panels. Width-fluctuation corrections were 
turned off when producing these results. 



B. Simulating the direct reaction of a surrogate 
experiment 

The first step of the reaction that takes place in a sur- 
rogate experiment is a direct reaction, such as a transfer 
or inelastic scattering reaction, that produces the rele- 
vant intermediate nucleus in a highly excited state. For 
the purpose of measuring the decay probabilities that are 
pertinent for the desired reaction, it is important that the 
intermediate nucleus first equilibrates into a compound- 
nuclear state. The relevant net result is therefore the 
distribution of JII states in the compound nucleus that 
is populated following the direct reaction. In principle, 



10 



300 



2T 200 




-I 1 1 — I I I I 



-I 1 1 — I I I I I I 1 1 T 



\ (a) 



Reference 



Decay model 1 

Decay model 2 

No wfc 



2.0- 

l-BK (b) 
1.6- 

1.4 

1.2 

1.0 



i — i i i 1 1 n | 1 i i i 1 1 i+f 



Decay model 1 
Decay model 2 
No wfc 



\ 



_i i i i i i ■ 



0.01 



0.1 



e n [MeV] 



FIG. 5: (Color online) Sensitivity of the 91 Zr(n, 7) 92 Zr cross 
section to variations of key parameters in the decay model, 
(a) Cross sections calculated using the three decay models of 
Table HQ plus one calculation using the Reference Decay model 
but turning off width-fluctuation corrections, (b) The ratio 
to the reference cross section. 



of a spin-dependent interaction, a distorted-wave Born 
approximation description of such a reaction yields only 



J" = 0+1-2+3- 



states. Distribution A, shown in 



Fig.[SJa), contains only such natural-parity states and we 
will use this distribution as the reference ("true") popu- 
lation of the compound nucleus 92 Zr* in a simulation of 
our benchmark surrogate reaction 92 Zr(a, a') 92 Zr*. For 
this reason we denote this distribution F A ' Tcl (J, II). Fur- 
thermore, the gamma-decay probability that is obtained 
when combining this JII distribution with our reference 
gamma branching ratio, G™ f (U, J, II), in Eq. ([5]) will be 
denoted P* cl (U). The decay probability P^ ci {U) thus 
corresponds to what would be measured in our simu- 
lated surrogate experiment, indicated as a solid curve 
in Fig.Ht a). 

In addition to F A ' rei ( J II) we use two more schematic 
distributions: F A+S {J, U) and F B (J, II). They are shown 
in Fig. EJb,c), respectively. The first one is almost iden- 
tical to F A ' Tei ( J, II), but contains a small random noise 
to simulate a minor error in the predicted distribution. 
Finally, Distribution B represents a significantly differ- 
ent population without any asymmetry between different 
parity states. The average angular momentum deposited 
in the compound nucleus, (J), is approximately 2.0 for all 
distributions. Therefore, these simulations correspond to 
a study of the sensitivity to deviations in the assumed JII 
distribution while keeping the average transferred angu- 
lar momentum approximately fixed. 



one would like to be able to describe this direct-reaction 
process, leading to an equilibrated compound nucleus, 
using an appropriate direct-reaction model. However, 
this is a non-trivial task since it requires a description of 
transfer and inelastic scattering reactions leading to un- 
bound states, as well as an understanding of the damp- 
ing of those states into equilibrated compound-nuclear 
states. Moreover, a variety of projectile-target combina- 
tions with a range of possible incident energies may be 
considered for producing the compound nucleus of inter- 
est. Different reaction mechanisms, regions of the nuclear 
chart, and projectile energies yield different compound- 
nuclear jn distributions and also provide different chal- 
lenges for a proper theoretical description. DWBA cal- 
culations relevant to the present work for inelastic alpha 
scattering to highly excited states in spherical nuclei arc 
currently under way and will be reported elsewhere. 

For the present purpose of simulating surrogate exper- 
iments we use simple schematic distributions. Further- 
more, we assume that the distributions are independent 
of the compound nucleus excitation energy in the en- 
ergy range just above the neutron separation energy that 
we are interested in. Since inelastic scattering reactions 
can be used to populate most compound nuclei that are 
relevant for studies of s-process branch points, we em- 
ploy distributions that exhibit the asymmetry between 
natural- and unnatural-parity states that is characteristic 
for those reactions with even-even targets. In the absence 



C. Three alternative analysis approaches 

In the following we will assume that Distribution A, 
F A ' re! ( J, n), is a reasonable representation of the JII 
distribution of a compound nucleus created in a surro- 
gate reaction (such as inelastic scattering off an even- 
even near-spherical nucleus). Also, we assume that the 
Reference Decay model with the parameters from Sec- 
tion IIVI represents a realistic description of the "true" 
decay of the compound nucleus, populated either through 
the neutron-induced or the surrogate reaction. The func- 
tion PJy Cl (U), calculated from Eq. ((5]) using the (energy- 
independent) formation probability f ( J, II) from Distri- 
bution A and the branching ratios G X {U, J, II) from the 
Reference Decay model, is then the quantity that is ob- 
served in our simulated Surrogate experiment. We now 
investigate various possibilities of extracting the desired 
"true" cross section cr^(e„), represented by the "Refer- 
ence" cross section in Fig. [SJ In order to test the different 
procedures, we introduce uncertainties in the modeling 
by using Decay model 1 or 2 (see Table ITT)) . This will il- 
lustrate the effect of having insufficient knowledge of the 
"true" decay of the compound nucleus. Furthermore, for 
the approaches that require the theoretical prediction of 
the JII population of the compound nucleus created in 
the surrogate reaction we make use of the schematic dis- 
tributions introduced in Sec. lVBl This will illustrate the 
effect of having insufficient knowledge of the jn distri- 
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FIG. 6: (Color online) Schematic surrogate populations used 
in this study: (a) F A ' Ic! (J, II); (b) F A+S (J,U); (c) F B (J,U). 
See text for further details. 



bution of the decaying nucleus. 

We discuss three different approaches that utilize decay 
data from a surrogate experiment to extract the desired 
low-energy (n, 7) cross section. The three approaches are 
labeled as follows: 

1. Weisskopf-Ewing approximation. 

2. Full modeling of the JII population following the 
surrogate reaction. 

3. Normalization of the decay model in the Weisskopf- 
Ewing region. 

Below we discuss each one of these approaches in some 
detail and draw some conclusions regarding their appli- 
cability. 



1. Weisskopf-Ewing approximation 



on the particular JII population of the intermediate nu- 
cleus due to the small number of open decay channels. 

Fig. [7] illustrates the inadequacy of the Weisskopf- 
Ewing approximation for this purpose. In this simu- 
lation the extracted (n, 7) cross section is obtained by 
simply multiplying the reference absorption cross sec- 
tion [cr^, (e n )] by the reference surrogate decay data, 
[P™ f (£/)]■ This procedure corresponds to the simulation 
labeled "WE" in Table HTT1 The extracted cross section 
is compared to the reference cross section. We find dis- 
agreement at the level of one order of magnitude, and also 
that the shape of the extracted cross section is wrong. 



TABLE III: Combinations of decay models (see Table HI)) and 
JII distributions (see Fig. |SJ used in the simulations. 



Simulation 



Decay model 



JII distribution 



Reference 
WE 
1 
2 
3 
4 



Reference 
Reference 

2 

2 

1 

1 



A,ref 

A+S 

B 
A+S 

B 



We briefly mention here the surrogate ratio method, 
which is based on the Weisskopf-Ewing approximation. 
This method has been used to extract (n, /) cross sec- 
tions in the actinide region [j| [f| . The goal of the ratio 
method is to determine the ratio of two cross sections of 
two "similar" compound- nuclear reactions. An indepen- 
dent determination of one of these cross sections then 
allows one to infer the other. The desired ratio of the 
two cross sections is obtained indirectly from surrogate 
measurements of the ratio of decay probabilities under 
the assumption that the Weisskopf-Ewing limit is valid. 
It was demonstrated in Ref. pj that the ratio approach 
may actually reduce the error that is usually associated 
with neglecting the JII dependence of the branching ra- 
tios. However, as we have just shown in Fig. the error 
encountered when applying the simple Weisskopf-Ewing 
approximation to low-energy (n, 7) reactions is about one 
order of magnitude larger than in the (n, /) case. In the 
latter case, the error rarely exceeds a factor of two even 
at small neutron energies (see in particular Fig. 10 of 
Ref. 0). Furthermore, for spherical or near-spherical 
targets, the JII dependence of the gamma-decay branch- 
ing ratios at low energies is very sensitive to the partic- 
ular level structure of the target, as shown in Sec. IIVI 
Thus, it is hard to reduce the effects of this sensitivity 
simply by measuring the ratio of decay probabilities for 
two neighboring isotopes. 



The results and discussion of the preceding sections im- 
ply that the convenient Weisskopf-Ewing approximation 
cannot be utilized in the surrogate approach when try- 
ing to extract low-energy (n, 7) cross sections for spher- 
ical and near-spherical targets. It was shown that the 
gamma-decay branching ratios depend very sensitively 



2. Full modeling of the JII population following the 
surrogate reaction 

For applications in the actinide region there have been 
attempts to take the effects of the JII population mis- 
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FIG. 7: (Color online) Extraction of (n,7) cross section from 
a surrogate experiment simulation using the Weisskopf-Ewing 
approximation; see simulation "WE" of Table Mil and discus- 
sion in text. 

match into account when estimating In, f) cross sec- 
tions from surrogate experiment data [ll], LL2|. In that 
work, the population of the compound nucleus follow- 
ing a (t,p) reaction was calculated in a distorted-wave 
Born Approximation approach. Furthermore, efforts are 
currently underway to develop a more advanced direct- 
reaction framework in which the surrogate reactions to 
unbound states, and their subsequent equilibration into 
compound-nuclear states, can be studied in greater de- 
tail. However, the application of the surrogate method 
to low-energy (n, 7) reactions is very challenging, in par- 
ticular when near-spherical targets are involved. We will 
show that a small uncertainty in the predicted JIT popu- 
lation can lead to a very large error in the extracted cross 
section. 

A surrogate analysis inspired by the diagram depicted 
in Eq. © is performed by first introducing a modeled 
gamma- decay probability P™ odol (£7), with both the JII 
distribution and the gamma-decay branching ratios ob- 
tained from initial modeling efforts 



pmodei^ = F raodol (c/, j, n)G^ odcl (c/, j n) 

J,n 

= ^]p™ odcl (c/, j,n). 



(12) 



This calculated quantity is then compared to the mea- 
sured gamma-decay probability. A fit to the experi- 
mental data is achieved by introducing a fitting function 
rj(U, J, IT) that relates the modeled decay probability with 
the measured one 



p^(u) = v(u, J, n)p™ dcl (u, j ii). 



(13) 



,/.n 



For simplicity we assume that the functions r](U, J, IT) 
are independent of J and II so that the fit corresponds 
to an energy-dependent normalization. In this case it is 
simple to evaluate the correction factor as 



n{U) 



P 7 cxp (E7) 

prnodel({7) ■ 



(14) 



and subsequently use it to extract the desired cross sec- 
tion 



<f act M = m J2 <b° dcl (£n, J, n) 



j.n 



x G™ dc \U, J,U)w ni (e n , J,n). 



(15) 



In the final step we have also introduced the theoretical 
width- fluctuation corrections, iu n7 (e n , J, II). The numer- 
ical calculations discussed below were calculated using 
the statistical-model code Stapre to calculate the var- 
ious factors in the above equation, including the width- 
fluctuation corrections, w ni {£ n , J II). 

We will now test this procedure by performing several 
simulations. We use the simulated surrogate decay data, 
P^ (U), as corresponding to the experimental data in the 
numerator of Eq. (fT4)) . As for the calculated quantity ap- 
pearing in the denominator, we introduce some variations 
in the modeled decay probabilities and the predicted JII 
population. The question is to what extent the analysis 
procedure will be able to correct for these variations. In 
the first simulation we use Decay model 2 to calculate 
G^ odel ([7, J,n). We know that this decay model overes- 
timates the cross section by -20 % (sec Table HI)- Fur- 
thermore, we assume that we have a very good, but not 
perfect, estimate of the JII population of the intermedi- 
ate nucleus. The distribution F modcl (U, J, II) only differs 
from the reference ("true") population F A ' tei ( J, II) by a 
small random noise. For this purpose we use the distri- 
bution F A+5 (J, II) shown in Fig.^b). The combination 
of Decay model 2 and the JII distribution F A+S (J, II) is 
denoted Simulation 1 in Table IIIII The procedure out- 
lined above results in the normalization function r](U), 
shown in Fig. Efb) plotted as a function of the equivalent 
neutron energy e n . Using this normalization function in 
Eq. (fTS"]) yields the extracted (n, 7) cross section shown 
as a dashed line in Fig. Oa). For this case the procedure 
works relatively well and the MACS at 30 keV for the 
extracted cross section is 67 mb. This number should be 
compared with the 74 mb that is the original prediction 
of Decay model 2, and the 62 mb that is the reference 
result that we were aiming for. 

In the second simulation we use a model for the JII 
population that is a poor representation of the "true" 
JII distribution. In this case, we find that the approach 
outlined above results in a very poor correction of the 
decay model, and consequently of the extracted cross 
section. For this simulation we use the combination of 
Decay model 2 and the JII distribution F B (J, II), which 
is denoted Simulation 2 in Table HITl The extracted (n, 7) 
cross section is shown as a dashed line in Fig. [5Jc), and 
the normalization function r)(U) is shown in Fig. [HJd). 
The extracted cross section is very different from the de- 
sired result below about 2.2 MeV. 

The above results imply that one may be able to use 
surrogate data to correct a theoretical decay model pro- 
vided one has sufficiently accurate information on the 
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FIG. 8: (Color online) Extraction of (n, 7) cross section from two different surrogate experiment simulations using a full 
modeling of the JII population: (a) Simulation 1, and (c) Simulation 2, see Table ITTT1 Panels (b) and (d) show the respective 
normalization functions r)(U). Note that U — s n + S„. 



JIT population of the compound nucleus that decays in 
the surrogate experiment. Obtaining a reliable predic- 
tion of the relevant JIT population is challenging, since 
it requires accurate direct-reaction calculations involving 
the nuclear continuum. Furthermore, the equilibration 
process that follows the production of a highly excited, 
intermediate nucleus in a surrogate reaction is not suffi- 
ciently well understood; possible decay mechanisms other 
than damping into the compound nucleus need to be ac- 
counted for if they are present. Nevertheless, it may 
be possible to obtain some experimental signatures of 
the JIT population of the decaying nucleus. While the 
high sensitivity of the gamma-decay branching ratios to 
the JIT population makes the extraction of the desired 
cross section from a surrogate experiment very difficult, 
it should also result in certain experimental obscrvablcs, 
such as the relative intensities of discrete gamma transi- 
tions, being useful to constrain calculated JIT distribu- 
tions. This remains to be investigated in more detail. 



3. Normalization of the decay model in the 
Weisskopf-Ewing region 

In the third approach we utilize the fact that surrogate 
experiments can, in principle, provide decay data for a 
very wide energy range. Thus it is possible to collect data 
from an energy region in which the Weisskopf-Ewing limit 



is approximately correct. For 92 Zr, we find from Fig. |H 
that this occurs at e„ approximately 3 MeV. A fit to the 
surrogate decay data in this region will be less sensitive 
to the predicted JIT population than at lower energies. In 
addition, the sensitivity studies for the calculated (n, 7) 
cross section, presented in Sec. IV Al showed that mod- 
eling errors in the s-wave average radiative width (r 7 ) , 
or in the level spacing Do, change primarily the mag- 
nitude of the calculated cross section, but do not affect 
the energy dependence of the cross section very much. 
Therefore, one may try to normalize the decay model in 
the Weisskopf-Ewing region, where the modeled quanti- 
ties are not very sensitive to the JIT population. The ex- 
tracted scaling factor is then used in the final calculation 
of the desired cross section at small energies. The normal- 
ization simply becomes an energy- and JIT-indepcndcnt 
factor, j]. 

We examine the outcome of this approach for three 
different simulations. As before, we employ the decay 
probability P^ cl (U) to represent the measured quantity 
in our simulated surrogate experiment. It is shown as a 
solid line in Fig.^a). The dashed line corresponds to the 
modeled decay probability obtained from Simulation 2, 
see Table IIII1 using the schematic distribution B (repre- 
senting a large error compared to the reference "true" 
distribution), combined with Decay model 2 (yielding 
gamma-decay branching ratios that are overestimated 
by ^20%). The dash-dotted line is the modeled decay 
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probability of Simulation 3, in which the schematic dis- 
tribution F A+S (J, IT) is used in combination with Decay 
model 1. In this case, the modeled JII population is very 
close to the reference distribution, but the decay model 
is one that produces gamma-decay branching ratios that 
are known to be a factor two too large. The dash-dash- 
dottcd line corresponds to Simulation 4 in which we still 
use Decay model 1, but the population of the interme- 
diate state is described by the schematic distribution B. 
From Fig. [Ha) we observe that the calculated gamma- 
decay probabilities differ by up to an order of magnitude 
in the most relevant energy regime. 

The next step in the analysis procedure is to con- 
struct the normalization functions T](U) as described in 
Eq. (fl4|) . The resulting functions are shown in Fig. [Hb). 
In this analysis approach, we focus on the high-energy 
region in which the normalization functions have become 
almost energy-independent. This is the region in which 
the Weisskopf-Ewing limit is approximately applicable. 
Still, we want to extract the normalization at an energy 
which is not too far away from the region of interest. In 
this particular case we use the energy e„ = 3 MeV, and 
for our three different simulations we obtain the normal- 
ization factors T) = 0.88, 0.50, and 0.54, respectively. The 
final step is simply to apply this constant renormaliza- 
tion to our modeling of the desired reaction in Eq. 
Again, we introduce width fluctuation corrections at this 
stage. Note that a different decay model was used in 
Simulation 2 as compared to Simulations 3 and 4, and 
that both of them differ from the reference ( "true" ) de- 
cay model. The final result of this procedure is shown in 
Fig. Wis) where the extracted cross sections are compared 
to the reference result represented by the solid line. The 
quality of the extracted cross sections is remarkable con- 
sidering the very different initial choices of decay models 
and JIT populations. The MACS at 30 keV for the three 
simulations are 65, 58, and 62 mb, respectively, to be 
compared to the 62 mb for the reference cross section. 

The underlying reason for the success of this approach 
is the direct proportionality between variations of the 
level density formula, or gamma-strength function, and 
the corresponding effect on the gamma-decay branching 
ratios as demonstrated, e.g., in Fig. [5] In the Weisskopf- 
Ewing limit this proportionality leads to a universal and 
energy-independent scaling of all JII components. This 
observation promises to be very useful for surrogate ex- 
periments in which data can be obtained for equivalent 
neutron energies at which the Weisskopf-Ewing limit is 
applicable. However, as we have seen in the 92 Zr ex- 
ample, see Fig. [9ja), the surrogate gamma-decay proba- 
bilities are quite small for energies above e n ~ 1 MeV, 
and consequently the measurement is challenging. The 
conclusion is that, given good-quality surrogate data in 
this energy region, the normalization approach outlined 
above offers an almost model- independent way of extract- 
ing the low-energy (n, 7) cross section for spherical and 
near-spherical nuclei. 



VI. CONCLUSIONS AND 
RECOMMENDATIONS 

We have examined the feasibility of using the surro- 
gate method to extract low-energy (n, 7) cross sections 
for spherical and near-spherical nuclei in the mass 90-100 
region. In particular, several Zr isotopes were studied 
and the 91 Zr(n, 7) 92 Zr reaction was used as a test case. 
Our study was performed by carrying out simulations 
to explore three different approaches to utilize surrogate 
reaction data. The sensitivity of the extracted cross sec- 
tions to uncertainties in the theoretical modeling was in- 
vestigated for the different approaches and their perfor- 
mance was assessed by comparing with a predetermined 
benchmark. 

One of the main results of the paper is the demonstra- 
tion of the large sensitivity of the gamma-decay branch- 
ing ratios, and consequently the surrogate-reaction 
gamma-decay probability, to the JIT population of the 
intermediate compound nucleus, and the effect that this 
sensitivity has on the extracted cross section. It was 
found that an approach in which the JIT population mis- 
match was not taken into account, gave an extracted 
cross section that deviated from the benchmark result 
by an order of magnitude. On the other hand, full mod- 
cling of the population of the compound nucleus is very 
difficult as it involves the description of direct reactions 
to continuum states and their subsequent equilibration 
into a compound-nuclear state. The third and final anal- 
ysis approach proposed in this paper promises to be the 
most viable one. It begins with a careful modeling of 
the decay, preferably using regional systematics to con- 
strain unknown nuclear properties. Surrogate decay data 
collected at slightly higher energies arc then utilized to 
normalize the modeled gamma-decay branching ratios so 
that the desired cross section can be extracted. 

It needs to be emphasized that these findings, and the 
sensitivities that were studied, apply to the analysis of 
our simulated surrogate experiment. In actual experi- 
ments, additional uncertainties are introduced, e.g., the 
identification of the final state, and the finite energy and 
angular resolutions. Furthermore, a prerequisite for our 
simulations of surrogate reactions has been the popula- 
tion of a fully equilibrated compound-nuclear intermedi- 
ate state. The question how this state is formed, through 
the population of highly-excited states in the direct sur- 
rogate reaction followed by multi-step equilibration pro- 
cesses, deserves to be studied in much greater detail. In 
particular, the probability for pre-equilibrium emission 
of particles should be investigated for different surrogate 
reactions. 

The theoretical modeling in this paper was per- 
formed using Hauser-Feshbach statistical-reaction the- 
ory. The approaches for normalizing the decay model 
using surrogate-reaction experimental data is based on 
the use of traditional decay models. Furthermore, these 
decay models were applied to describe the decay from 
all relevant JIT states in the entire energy range of in- 
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FIG. 9: (Color online) (a) Gamma-decay probabilities for different simulations. The solid line corresponds to the reference model 
whereas the other three are based on different combinations of the assumed JII population and decay model (sec Tabic [TlTTl . 
(b) The normalization functions r/(U). Their value at e„ = 3 MeV is used in the final cross section calculation. Note that 
U — En + Sn- (c) The extracted cross sections, for the respective simulations, obtained by employing the analysis approach 
outlined in Sec. IVC3I The solid line is the reference cross section. 



terest. This implies that the arguments presented here 
are not applicable to regions of nuclei that are too far 
away from the valley of stability, where the use of estab- 
lished level-density formulae and parameterized gamma- 
strength functions has to be questioned. The same obser- 
vation applies to the optical potentials that were used in 
this work to calculate particle transmission coefficients. 
The normalization procedures outlined in Sec. IV Cl rests 
on the accuracy of the optical model that is used to 
compute the absorption cross section. Finally, we stress 
that width fluctuation corrections were always applied 
in the final modeling part of the analysis. Those correc- 
tions cannot be obtained from the surrogate experimental 
data. 

With the insights gained from this study we can give 
a few general recommendations for future surrogate ex- 
periments: (1) Good particle identification is very im- 
portant since one needs the absolute normalization to 
get the decay probabilities with good accuracy. We note 
that the gamma-decay probability is expected to decrease 
dramatically in magnitude and one must be able to mea- 
sure them up to a few MeV above the neutron separation 
energy. (2) Consequently one also needs to be able to ef- 
ficiently identify the desired decay channel. For gamma 
decay this might require the existence of a very strong 
collector 2 + — > + transition to tag the occurrence of 



a gamma cascade. This requirement would limit the 
method to cases where the intermediate nucleus is an 
even-even isotope. However, other options for tagging 
the gamma cascade are being considered (37}- (3) The 
energy resolution is less of an issue if the approach of 
renormalizing the decay model at high energies is used. 
In this procedure the energy dependence is obtained from 
the statistical-reaction modeling. Still, an energy resolu- 
tion that is better than 100 keV should be desirable. (4) 
Additional experimental information, such as relative in- 
tensities of the gammas de-exciting low-lying levels of 
different spin and parity, can help to gain insight into 
the JFl population of the compound nucleus. However, 
one should remember that the observed gamma intensi- 
ties depend on the properties of the gamma cascade that 
proceeds through the quasi-continuum as well as on the 
initial JII distribution, and so it is important that the 
properties of this cascade be accurately modeled. (5) 
Measuring the decay at different emission angles of the 
ejectile in the initial direct reaction should provide an ex- 
perimental handle to vary the JII population. However, 
it remains to be seen from direct-reaction modeling for 
each particular case how large this change can actually 
be. (6) Finally, we recommend that benchmark experi- 
ments, in which the surrogate data is used to extract a 
known cross section, are carried out since this will pro- 
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vide very valuable insights into the issues discussed in 
this paper. 

Finally, we note that applications involving spherical 
nuclei probably constitute one of the most challenging 
applications of the surrogate method. This statement ap- 
plies, in particular, to isotopes near magic shells such as 
the Zr isotopes studied in the present work. In addition, 
low-energy (n, 7) reactions provide the most difficult re- 
action with regards to the angular momentum mismatch. 
In contrast, for deformed nuclei there is usually a larger 
number of open decay channels and that reduces the sen- 
sitivity of the gamma-decay branching ratios to the initial 
JII population. In particular, the Weisskopf-Ewing limit 
might be reached already at relatively low excitation en- 



ergies above the neutron separation energy 
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